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Abstract 

In a recent paper |lj we found strong evidence from simulations that the Ising antiferromagnet on 
"thin" random graphs - Feynman diagrams - displayed a mean-field spin glass transition. The intrinsic 
interest of considering such random graphs is that they give mean field theory results without long 
range interactions or the drawbacks, arising from boundary problems, of the Bethe lattice. In this 
paper we reprise the saddle point calculations for the Ising and Potts ferromagnet, antiferromagnet 
and spin glass on Feynman diagrams. We use standard results from bifurcation theory that enable us 
to treat an arbitrary number of replicas and any quenched bond distribution. We note the agreement 
between the ferromagnetic and spin glass transition temperatures thus calculated and those derived 
by analogy with the Bethe lattice or in previous replica calculations. 

We then investigate numerically spin glasses with a ±J bond distribution for the Ising and 
Q = 3, 4, 10, 50 state Potts models, paying particular attention to the independence of the spin 
glass transition from the fraction of positive and negative bonds in the Ising case and the qualitative 
form of the overlap distribution P(q) for all of the models. The parallels with infinite range spin glass 
models in both the analytical calculations and simulations are pointed out. 
Submitted to Nucl. Phys. B. [FS] 



1 Introduction and Analytical Calculations 



The analytical investigation of spin glasses on random graphs of various sorts has a long and honourable 
history |^|, [| , though there has been little in the way of numerical simulations Random graphs with a 
fixed or fixed average connectivity have a locally tree like structure, which means that loops in the graph 
are predominantly large, so Bethe-lattice-like Q (ie mean field) critical behaviour is expected for spin 
models on such lattices. Given this, the analytical solution for a spin model or, in particular, a spin glass 
on a Bcthe lattice Q can be translated across to the appropriate fixed connectivity random lattice. 
Alternatively, a replica calculation can be carried out directly in some cases for spin glasses on various 
sorts of random lattices. For the case of Ising spins on a fixed connectivity random lattice we arrive at 
the following prediction for the ferromagnetic transition temperature, if one exists 

P(J) tanh(/3 FM J) dJ = (1) 

where P{J) is the quenched probability distribution of bonds in the model, z is the connectivity and 
Pfai is the (inverse) transition temperature. A similar relation is predicted for the spin glass transition 
temperature (3sg 

[■CO i 

/ P(J) tanh 2 (/3 SG J) dJ = ——. (2) 

From equs.(|]J|) we would expect to find a ferromagnetic transition when P( J) = p 6(J— 1)+(1— p) S( J+l) 
at 

(2p - 1) tanh(/W) = ( 3 ) 

or a spin glass transition at 

tanh 2 03 SG ) = (4) 

depending on which critical temperature was lower. The global order parameter that is necessary to 
describe the system in the spin glass case [|[ [| appears as the fourier transform of the local field distri- 
bution. 

A rather different way of looking at the problem of spin models on random graphs was put forward in 
||] , where it was observed that the requisite ensemble of random graphs could be generated by considering 
the Feynman diagram expansion for the partition function of the model. For an Ising ferromagnet with 
Hamiltonian 

H = pJ2 a * a 3> ( 5 ) 

<ij> 

where the sum is over nearest neighbours on three-regular random graphs (ie 4> 3 Feynman diagrams), the 
partition function is given by 

where N n is the number of undecorated graphs with 2n points, 

A\ 2 " ( 6n-1)!! 
n ~U/ (2n)!! [) 

K is defined by 

K ab = \ V V £ ) (8) 

1 For some recent simulations, sec jij. 
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and the action itself is 



s = lj2<i>aK-b<Pb-l(<t>l + <£-)■ (9) 

a, b 



where the sum runs over ± indices. The coupling in the above is g = exp(2/3J) where J = 1 for the 
ferromagnet and the </>+ field can be thought of as representing "up" spins with the cf>- field representing 
"down" spins. An ensemble of z-regular random graphs would simply require replacing the </> 3 terms with 
(j) z and a fixed average connectivity could also be implemented with the appropriate choice of potential. 

This approach was inspired by the considerable amount of work that has been done in recent years on 
N x N matrix ^ versions of such integrals which generate "fat" or ribbon graphs graphs with sufficient 
structure to carry out a topological expansion |lC| ] because of the matrix index structure. The natural 
interpretation of such fat graphs as the duals of triangulations, quadrangulations etc. of surfaces has 
led to much interesting work in string theory and particle physics The partition function here 

is a poor, "thin" (no indices, so no ribbons), scalar cousin of these, lacking the structure to give a 
surface interpretation to the graph. Such scalar integrals have been used in the past to extract the 
large n behaviour of various field theories [ ft2| again essentially as a means of generating the appropriate 
Feynman diagrams, so a lot is known about handling their quirks. 

For the Ising ferromagnet on three-regular (</> 3 ) graphs, solving the saddle point equations at large 
n shows that the critical behaviour appears as an exchange of dominant saddle point solutions to the 
saddle point equations 

4>~ = Vg0- + ^l (io) 

at g = exp(2/3pM) = 3. The high and low temperature solutions respectively are 

_ V9 



.9+1 



^ fi±,^EF] ,ii, 



2(5-1) V V .9 + 1 

which give a low temperature magnetized phase that can be detected by a non-zero magnetization order 
parameter, which is defined in terms of the fields cf> + , 4>- as 

M = \+ \~ . 12 

The critical exponents for the transition can also be calculated in this formalism and, as expected, are 
mean field. In general a mean field transition appears at 

exp(2l3 FM ) = z/(z - 2) (13) 

on (j) z graphs, which is exactly the value predicted by equ(||) in the standard approaches. 
To consider a non-trivial bond distribution it suffices to make the substitution ||] 

/oo 

Kab P{J) dJ. (14) 
-oo 

in the saddle point equations. If we take the distribution P(J) = p 6(J — 1) + (1 — p) S(J + 1), this has 
the effect of replacing g by 

P9+1-P (15) 

P+(1-P)9 (5) 



"N, the size of the matrix is not to be confused with n, the number of vertices in the graph! 
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On (j) z graphs we find that this change in the coupling shifts the ferromagnetic transition point to 



(2p-l)(z-l) + l 
6XP(2/W) = (2p-l)(z-l)-l (16) 

which is again identical to the value found from equ.(||). 

This is rather surprising as a direct substitution of the weighted propagator coming from equ. ( |l4| ) into 
the saddle point equations might be expected to correspond an annealed distribution of bonds because 
we are calculating with a finite number (one) of replicas. The calculations leading to equ.(||) and equ.(Q), 
however, have taken quenched bond distributions by calculating with k replicas and then taking the limit 
k -> 0. 



2 The Saddle Point Equations for Spin Glasses 

There are various possibilities for addressing spin glass order in the Feynman diagram approach. In 
the entropy per spin was calculated for the Ising cmfi-ferromagnet on </> 3 graphs and it was found to 
become negative for sufficiently negative /?, which is often indicative of a spin glass transition. Similarly 
it was found that the factorized solution (which exists in the ferromagnetic and anti-ferromagnetic cases) 
broke down for higher moments Z k of the partition function, which is again indicative of a spin glass 
transition. The temperature at which this happened appeared to be converging to a finite value as the 
moments increased, unlike the random energy model where a similar calculation gives a temperature that 
diverges as \fk |p^| . 

We now look at the calculation of Z 2 in a little more detail. We can put two Ising models on each 

Feynman diagram by taking four fields , </>h , 4> v > 4> i where the double subscript now covers both 

replicas, and enlarging the inverse propagator to 
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(17) 



If we solve explicitly the saddle point equations for the two replica system with P(J) = pS(J — 1) + 
(1 — p)S(J + 1) we find a p-independent solution appearing at g = 5.8284.. and a mirror solution at the 
corresponding antiferromagnetic value given by the inverse of this, g = 0.1716.... The second solution, 
g = 0.1716..., was also observed in M as the point where the factorized solution broke down for the 
antiferromagnet with k = 2. Remarkably, as we are looking at a finite number (two) of replicas here, 
g = 5.8284, g = 0.1716 are precisely the values given by eqn.(^) for /3sg, where the limit k — > has been 
taken to obtain the spin glass transition temperature. Thus, just as for the ferromagnetic transition, 
the quenched k — > results are appearing already at k ^ 0. Using numerical routines to investigate the 
structure of the saddle point equations for more Ising replicas reveals that the solutions at 0fm and (3sg 
still appear at higher k as bifurcations from the symmetric high temperature solution. 

This is not quite the whole story. As already noted in || for k = 3,4... we hit a different sort of 
critical behaviour which takes us off the symmetric solution branch to a replica symmetric solution before 
reaching Psg- At these points a first order (jump in q) transition occurs. This behaviour is identical 
to that of the Ising replica magnet [jl4j , which is effectively the SK model at a finite number of replicas 
for a suitable choice of parameters. In this the k — 2 transition is second order and the transition 
temperature coincides with the k = value, whereas the k > 3 transitions are first order and occur at 
higher temperatures. 

We can attempt to analyse the saddle point equations for any number k of Ising replicas in the case of 
purely ferromagnetic or antiferromagnetic couplings using bifurcation theory methods because the tensor 
product structure of the propagator K a b is preserved in these cases which facilitates the calculation. 

3 This is just a more formal statement of the observation in |J that a solution given by the product of k = 1 solutions 
existed in these cases. 



3 



We will not, however, see first order transitions such as those discussed above because the bifurcation 
structure only tells us about continuous transitions from the symmetric high temperature solution. 

If we denote the 2x2 propagator in equ.(7) as K{g) for brevity, the saddle point equations for the k 
replica case may now be written schematically as 

5^^-®^ = (18) 
ocp 

where (f>, <j) 2 are now 2 k dimensional vectors, the tensor product is taken over k copies of K and the action 
for k replicas is denoted by S}.. We have abused the notation slightly in equ.(18) as the left hand side is 
really ® k K(g)dS k /d4> because we have multiplied by ® k K(g) to get rid of the inverse factor coming from 
the quadratic term in S k , which simplifies the analysis of the saddle point equations. We will continue 
the abuse in the discussion below. 

For the antiferromagnetic case g is simply replaced by 1/g throughout. We restrict ourselves to 
investigating the symmetric branch (and hence ignore the first order solutions that appear for k > 2). 
A symmetric, high temperature solution can be constructed by taking products of the k — 1 symmetric 
solution in equ.(10) (t,t) where 

t=^-. (19) 
.9 + 1 

Standard theory [jl5| then shows that a bifurcation from the symmetric high temperature solution (ic 
phase transition) is expected when the Hessian det (^d 2 Sk/d$ 2 ^j is equal to zero. More explicitly 

det ( ^-^] = det(2t k K k - 1) (20) 



where we have denoted the tensor product of k copies of K by K k - If we absorb the factors of t into K 
and denote K = (yfg + 1/^/g) -1 K we thus have 

det ( ^-^] = det(2K k - 1) = (21) 



as our bifurcation equation. The extension to z-regular (0 Z ) graphs is trivial. For these the k = 1 high 
temperature solution is (i, t) where 

f. (^)'"'" W 

and the Hessian becomes 

det (H^) = det((z - lW- 2 ) k K k - 1) (23) 

so the net effect is simply to replace 2 by z — 1 in the bifurcation equation. 

The matrix K has two eigenvalues 1 and (g — l)/(g + 1) = tanh(/3). From the algebraic structure of 
K k we can readily calculate the eigenvalues A m = tanh(/3) m , m = 0, 1, 2, . . . k where each eigenvalue has 
multiplicity ( ) and express (taking z — 3) 

det = fl (2tanh(/?r - 1)® = 0. (24) 

The roots of the bifurcation equation are thus 

tanh(/3) = 2- 1/m m = 1, 2, . . . fe (25) 
reproducing equ.(3) for p = 1 and z = 3. In terms of /3 = (1/2) log(g) we can write this as 

1, /2 1 /' m + l 



lm = 5 log( 2l/m _ 1 ) m=l,2...A (26) 
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where we have retained the label m on (3. The above analysis holds for the ferromagnet where g = 
exp(2/3) > 1. In the case of a pure antiferromagnet g = exp(2/3) < 1 and tanh(/3) is negative, so we can 
only find roots of the bifurcation equation for m even. 

1 /2 1 A»_i 



^=2 l0 H2TMTlJ - = 2,4. ..fc (27) 

We have thus recovered both (3fm = Pi (g = 3) for the ferromagnet and (3sG = @2 (<? = 0.1716) for the 
antiferromagnet in this approach. 

A generic feature of the above solutions is immediately obvious - once a bifurcation appears at some 
k it is always present for a greater number of replicas. This suggests a possible resolution of the puzzle of 
why the quenched transition temperatures are appearing already at finite k: The appropriate transition 
is encountered at the same temperature for all k, so the analytical continuation k — > is trivial. In the 
spin glass case the first order transitions for k > 3 are the fly in the ointment, as these are encountered 
before the putative continuous transition at (3sg- Looking explicitly at the replica symmetric solution for 
k = 3 Q shows that the the replica symmetric banch appears in a first order transition at g ~ 0.19 . . . 
(in the antiferromagnet) before the symmetric branch takes over again below this. 

It is rather amusing that we have encountered a "c=l (2 Ising model) barrier" in this context. Calcu- 
lations are still possible here for k > 2, but there is certainly a change of behaviour. In the case of spin 
models on annealed ensembles of fat graphs analytical calculations have struggled to get beyond c = 1 , 
whereas simulations have so far failed to see much difference between c < 1 and c > 1. 

Although the tensor product structure of is lost when p ^ 0,1 and the saddle point equation 
becomes 

^ = $- (p ® k K(g) + (1 - p) ® k K(l/g)) $ 2 = (28) 

the matrices K(g),K(l/g) differ only in the sign of the eigenvalue tanh(/3), so it is still possible to derive 
an explicit expression for the Hessian, which in this case is 

det f5r) = II [2tanh(/3r( P+ (-ir(l-p))-l] a) - (29) 
V or / m=0 

We thus reproduce the values of (3m and Psg calculated from equs.(3,4) for general p as the first two 
bifurcation points. On cj> z graphs we replace 2—*z—\, which is still in accordance with equs.(3,4). As all 
even bifurcation points are independent of p, we predict the p independence of the spin glass transition 
point, as in the more standard approaches. 

Equs.(3,4) are not restricted to distributions of the form P{J) — pS(J — 1) + (1 — p)S(J + 1), but we 
can extend the approach here to accommodate this by making the observation that K(g\) will commute 
with another if (52) for any 31,32- Thinking of a given distribution (a gaussian, say) as a weighted sum 
of delta functions over if's with different values of g, we can still diagonalize all of the terms in the sum 
simultaneously and arrive at an expression of the form 



M5r) = n 



m— 



P(J) tanh(/3J) m dJ - 1 



lit) 



= 0. (30) 



for the bifurcation equation. 

Although the continuous transition temperatures (3 m are independent of the number of replicas, the 
multiplicity of solutions (ie low temperature phases) that bifurcate at a given point are not. There are 
fc(fc— l)/2 solutions bifurcating at the spin glass transition temperature /32, so just as in the infinite range 
model there is a change in the nature of the solutions to the saddle point equation when k < 1 and the 
number of solutions is formally negative. The classification of the spin glass solutions in the current case 
is thus clearly closely related to the parametrization of the matrix Q a b that appears in the saddle point 
equations of the infinite range model |lq] . 
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It is also possible to consider Potts spins on a Bethe lattice |?], P . We take Hamiltonian in this case 
to be 

= 2/3 £ ^ ( 31 ) 

<ij> 

where the spins <7j can now take on Q values 0. The critical temperatures for ferromagnetic and spin 
glass ordering calculated in @] give, in the style of equs.([l],||), 

/ exp(2/W)-l \ dJ =^_ (32) 
, V ; \exp(2p FM J) + Q-l) z-1 V ' 

and 

/ exp(2fe G J) - 1 V 1 

F{J) \cMWsgJ) + Q-i) z l ^ 

for Q state Potts models. If we consider the ferromagnet, P(J) = 5(J — 1), we find that exp(2/?i?M) = 
(z + Q — 2)/(z — 2). There is, however, a parallel here with the Ising spin glass calculations. It was noted 
in || that this temperature is actually a spinodal point, for Q > 2 the Potts ferromagnet undergoes a 
first order transition before reaching this point. Arguing, as in the Ising case, that a random graph looks 
locally like the Bethe lattice we would expect the various critical temperatures and spinodal points still 
to apply on <fi z graphs. 

It is also possible to solve the saddle point equations explicitly to avoid relying on this analogy. For 
the 3-state Potts model with action 

S = +$ + €) - C&afa + + Mc) ~ + 4 + 4>% (34) 

this gives high and low temperature solutions^] 

4>a,b,c = 1 - 2c ; 



1 + Vl - 4c- 4c 2 



1 + 2c - Vl - 4c - 4c 2 



(35) 



where c = l/(g + 1). This gives a changeover in behaviour at c = 1/5, ie g = 4, which agrees with the 



value from equ.(32) and thus picks out the spinodal point. Similarly, solving the 4-state Potts model 



equations gives the solutions 

(t>a,b,c,d 



1 - 3c 



l + 3c- Vl-6c-3c 2 



,b,c 



l-c + Vl-6c-3c 2 
<f>d = g ( 36 ) 

where c = 1/(^ + 2) here. The changeover here is at g = 5, again apparently at the spinodal point. It 
would be interesting to carry out a simulation to check the nature of the transition in these ferromagnetic 
models. 

Analysis of the Potts spin glass saddle point equations follows a similar tack to those of the Ising 
model. The caveats about missing potential first order transitions at finite k still, of course, apply. If we 

4 We have adopted a different normalization for the Potts Hamiltonian from [fl, taking a factor of 2 rather than Q in 
front, for consistency with \ 



5 We gave these, and those for the 4-state model, incorrectly in [0. 
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take the three state Potts model as an example K is replaced by the 3x3 matrix L(g) 

Lab = 

and the high temperature solution on <f> 3 graphs is (t, t, t) where 




(37) 



(38) 



The eigenvalues of L(g) = tL(g) are 1, A, A, where A = (g — l)/(.g + 2). It is still true that L(gi), L(g2) 
will commute for any gi, gi so we can put this to good use again to derive the Hessian for general P(J) 



det 



d 2 S k 
d$ 2 



n 



v ' Vexp(2/3J) + 2 



UK 



(39) 



The multiplicity of solutions has changed because of the possibility of taking one of two A from each 
replica. The results for /3fm and /3sg from the Bethe lattice are thus faithfully reproduced as the first 
two bifurcation points. 

Generalizing to the Q state Potts model on <f> z graphs, we find 



t = 



9 



g + Q-l 



l/(z-2) 



(40) 



and the eigenvalues of the Q x Q matrix L will be 1 and Q 
gives the Hessian 



1 A's, where A= (g-l)/(g + Q- 1). This 



det 



d 2 S k 



n 

m=0 



(Z-1) P(J) 



( exp(2/3J) - 1 
Vexp(2/3J) +Q- 



dJ -1 



(™)(Q- 



(41) 



which the reader will be pleased to hear is the last elaboration of the bifurcation equation that we consider. 
The Q = 2 case reproduces the Ising results derived earlier. 

Given the accessibility of mean-field results on random graphs, especially as one is not forced to 
take infinite range interactions or deal with the boundary difficulties of the Bethe lattice, the lack of 
simulational effort is rather surprising. As the new viewpoint offered by regarding the random graphs 
as Feynman diagrams offers the possibility of attacking various questions, such as the nature of replica 
symmetry breaking on random graphs, from a different angle, we thought it important to continue the 
simulations of jlj] for the Ising and Q state Potts spin glasses on cj) 3 graphs to check carefully the agreement 
with the known results for quantities such as the transition temperatures and critical exponents and to 
gain a thorough understanding of the numerics before going on to investigate less well understood aspects. 
Our numerical methods are largely similar to those of Q, so we describe them only briefly in the next 
section and concentrate on the results of the simulations. 



3 Simulations: Generalities and Ising Model 

In all of the simulations we report here we generated around 100 different Feynman diagrams at each 
of the sizes we simulated, 100, 250, 500 and 1000 vertices. The quenched bond distribution on each 
graph was P{J) = p S(J - 1) + (1 - p) 5{J + 1), with p taken to be 0.5 (the "±J" distribution) 
unless specified otherwise. Various massively parallel processors were used to run the simulations of the 
different graphs simultaneously and perform the quenched averaging over the graphs in situ. Again, 

6 This definition of L. which is the simplest from the point of view of the saddle point equations, involves a rescaling 
with respect to equ.(BJ) which accounts for the different between t and the high temperature solution in equ.(p5|). 
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in the interests of simplicity and reliability (and as in |D) we used only a metropolis algorithm with 
simulated annealing rather than some of the more advanced algorithms developed recently. At each 
(3 value simulated we carried out 500,000 simulated annealing sweeps, followed by 500,000 production 
sweeps with a measurement every tenth sweep. Each sweep consisted of a complete metropolis update 
of the lattice. Our strategy for extracting the critical exponents will be to use Binder's cumulant for the 
overlap, defined below, to extract an estimate for (3$g and the combination vd and then to look at the 
finite size scaling of other quantities to extract further exponents. 

In the Ising model the spin glass transition temperature is found in simulations by putting two Ising 
replicas on each graph with spins er^, ti to measure the overlap, which is the order parameter for the spin 
glass transition 

g=-Vcr l r l . (42) 
n £ — ' 

i 

The Binder's cumulant for the overlap is then defined in an analogous fashion to the Binder's cumulant 
for the magnetization in a ferromagnetic transition 



j(i 4 ) 
K? 2 



Us 9 = (43) 



where ( ) denote thermal averages and [ ] disorder averages. The plots of U sg for differently sized graphs 
are expected to cross at Psg- In the course of all the simulations the overlap distribution 

P w («)= [<%-l/n5>i7-<)>] (44) 

was also histogrammed. A non-trivial P(q) is a strong, but not infallible, signal for a spin glass phase. 
It is perhaps worth remarking that the order of averages in eqn.(^3|) is not absolutely obvious a priori 
- one might have considered [ < q 4 > / < q 2 > 2 ], for instance. Regarding < q n > as moments of 
the distribution P(q), the choice in the eqn. (f43|) would seem most appropriate. The alternative overall 
average has been considered in quantum spin glass simulations |TjJ on the heuristic grounds that it gives 
better scaling behaviour. We have carried out the scaling analysis of our simulations described below 
with both definitions for the Ising spin glass and found essentially identical results. 

We expect the overlap to have the finite size scaling form q ~ n -Pl vd at any spin glass transition and 
the spin glass susceptibility, defined as 



\Y>^\ ( 45 ) 

or, alternatively, 



Xs 3 n 



Xsg =n J q 2 P n (q)dq. (46) 

to diverge as \sg — n 1 l vd . Other critical exponents such as that for the specific heat C ~ B + Con a / vd 
are defined in the standard manner and relations such as a = 2 — vd still hold good. 

The mean field critical exponents are shown in the table below for an Ising (Q = 2) spin glass 

Table 1: Spin Glass exponents for Ising model 
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We have written vd in the table rather than v as we cannot disentangle v and d in the infinite dimensional 
random graph case at hand. 

We now commence the analysis of the results of the simulations proper for the Ising model. The 
cumulant for the overlap is plotted using both ways of carrying out the averages in Fig.l and Fig. 2, from 
which it is clear that the resulting crossover points are, as claimed, very similar. Fitting the crossover 
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point from either figure gives Psg — 0.88(1), although the overall average that gives Fig. 2 is slightly 
noisier. This value is in good agreement with the theoretical predictions. An additional consequence of 
equ.(0) is that (3sg should be independent ofp for P(J) — pS(J— 1) + (1— p)6(J+l). In Q we considered 
the Ising antiferromagnet, which has p = 0, and found (3sg — —0.94(2) by direct observation of cumulant 
crossing, but we also found (3sg = —0.88(2) from extrapolation of the specific heat peak. These results 
are thus only marginally compatible with the p independence. As a further check we simulated a complete 
set of graphs with p = 0.3 and found that the cumulants were identical within the errors to the p = 0.5 
values, this giving the same Psg- The slight discrepancy at p — may well be due to the poorer statistics 
we have for these runs. 

Further analysis of the cumulants allows us to extract vd by looking at the scaling of the maximum 
slope of the cumulant with the number of vertices n. 

f du \ — 

max — ~ ri"<J (47) 



dp 

We find that vd = 2.9(1) from the data in Fig.l and vd — 3.0(2) from that in Fig. 2. The simulations in 
(!]] found the compatible value of vd — 2.8(2) for p — and all are in agreement with the mean-field value 
of a = — 1 for the specific heat critical exponent, deduced from a = 2 — vd. A cusp in the specific heat 
C, rather than a divergence, is indicated by this negative value of a in the finite size scaling relation 

C ~ B + C Q n a ' ud (48) 



The cusp can clearly be seen in Fig. 3. A direct fit to the scaling of the specific heat as in eqn.(|48|) is 
not particularly instructive, as the extra adjustable constant B allows for a very good fit to the value 
a/vd = —1/3, but it is at least consistent with the value emerging from the cumulant analysis. 

Turning now to the expected scaling of the overlap, we find respectable agreement with the finite size 
scaling relation 

q ~ n-Pl vd (49) 

with (3/vd calculated to be 0.30(1) by extrapolating to the pseudocritical value of (isc = 0.88. Having 
considered the first moment of P(q) with q, we now move on to what is effectively the second moment of 
the P(q) with Xsg- In Fig. 4 we see the expected divergence, and a fit to the finite size scaling relation 

Xs g * ntl" d (50) 

gives j/vd — 0.34(2), in agreement with the mean field value of 7 = 1. No divergence is expected in the 
linear susceptibility 

dM (n\ 
XM = m (51) 

and, as is clear from Fig. 5, none is seen. We have not measured the response of the model to an external 
field, nor attempted to fit the correlation functions at the critical temperature, so we have no estimates 
for 5 and rj. However, it is clear from the above results that the model is giving mean-field like behaviour 
for the spin glass transition that is observed. 

The results for the specific heat and the susceptibilities reported here for the ± J spin glass are essen- 
tially identical to those in jlj for the Ising antiferromagnet, which can be taken as a further confirmation 
of the p independence of the spin glass transition. Looking directly at the histograms of P{q) themselves 
in the spin glass phase for p = (the antiferromagnet) and p = 0.3, 0.5 (the models simulated here) does 
show some differences as can be seen in Fig. 6. However, the determination of the critical point depends 
on the scaling of the moments of P(q) around the transition point at (3 = 0.8814... and it is here we 
expect to find p independence. This is indeed what we see at j3 = 0.9 in Fig. 7 - the histograms are 
essentially identical for p — 0.3, 0.5. There is a slight discrepancy still for p = 0, which is a reflection 
of the less convincing scaling analysis in jij that determined the critical point in the case of the pure 
antiferromagnet. We have shown only the P{q) for graphs with 500 points in the figures for clarity, the 
results for other graph sizes are identical. 
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4 Simulations: Potts Models 



For Q-state Potts models the overlap may be defined as 

1 " 

q = -Y,(Q8* i ,T i -i) (52) 

n L — ' 

i=l 

which is arranged to be zero in the uncorrelated case. The cumulants and P(q) may then be dehned in 
an analogous fashion to the Ising model. There is, however, no longer the q — > —q symmetry in P(q) that 
is present in the Ising spin glass for Q > 3 state Potts models. In all of the simulations reported here 
we took P(J) = 1/2(5(J - 1) + S( J + 1)). We simulated lattices of size 100,250,500 and 1000 for the 
Q = 3, 4 state models, but only the three smaller lattice sizes for Q = 10, 50. The statistics and annealing 
schedules were identical to the Ising model simulations. We are principally interested in observing the 
qualitative features of the various Potts models, so our analysis concentrates less on serious finite size 
scaling to extract the exponents than the Ising results. 

Mean field theory (in the infinite range model) suggests that Q > 3 Potts glasses are rather different 
in behaviour from the Ising spin glass as there is no longer continuous replica symmetry breaking. For 
Q = 3, 4 there are two consecutive transition temperatures PsG and PsG2 to different glass phases and for 
Q > 4 the first has a discontinuity in the overlap q. For f3sG2 > (3 > (3sg P(q) consists of a delta function 
at q = and another delta function at finite q with no continuous features whereas for j3 > PsG2 the delta 
function at finite q splits in two with a continuous distribution between. The values of the transition 
temperature given by equ. (33,41) are for ftsc- Some features of the solution differ on the Bethe lattice 
, in particular an Ising- like solution appears to exist for coordination number three with a ± J bond 
distribution, and the <f> 3 random graph model might be expected to behave in this fashion too. 

If we take the 3 and 4 state Potts models first, we can look at the cumulant crossing in a similar 
manner to the Ising model to attempt to pinpoint the phase transition. It is possible to extract crossing 
points at /3 = 1-5(1) for Q = 3 and f3 — 2.0(1) for Q = 4, which are in agreement with the values 
calculated from equ.(|33|) shown below in Table. 2. Our data points for the Q = 10,50 state models are 
rather sparser but the cumulant crossing appears to become more clear cut with increasing Q, which 
offsets this. 



Q 


2 


3 


4 


10 


50 


0SG 


0.8814 


1.5824 


2.1226 


3.6914 


6.1943 


measured 


0.88(1) 


1.5(1) 


2.0(1) 


3.6(4) 


6.1(1) 



Table 2: (3sg f° r various Q state Potts models on </> 3 graphs. 

The quoted errors in the table are the most conservative choice - the largest difference between the 
crossing points for the various lattice sizes - apart from Q — 3 where the N — 100 data was dropped, as 
it failed to cross with the other lattice sizes. 

We next look at the qualitative features of the low temperature phase in the various P(q) to see if the 
mean field theory expectations are borne out. With our definitions, the histograms of P(q) range from 
— 1 to Q — 1 for a Q state Potts model with the origin at zero and roughly half of the probability density 
in the range 0. . . (Q — 1), so we have plotted only ... (Q — 1), normalized to the range ... 1 for clarity 
in all our figures [] In all of the models the high-temperature P(q) look like skewed gaussians centred 
on the origin, whose sharpness increases with graph size indicating delta-function-like behaviour in the 
continuum limit. In Fig. 8 we have plotted P(q) for the three state Potts model for (3 — 2.2 (> Psg)- 
This is clearly different from the Ising model: there is a strong peak at the origin that is increasing with 
graph size as well as smaller bump at around q = 0.7 which is also increasing with graph size and which 
might tentatively be identified with another peak in the continuum limit. In Q simulations of the Potts 

7 It is also possible to extract a positive scalar order parameter for the spin glass phase by going to a simplex representation 
for the Potts spins and using the radial component of the overlap matrix, but we stick with the canonical definition here. 
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antiferromagnet, where no spin glass phase is expected, saw no signs of the structure at q = 0.7, with 
P(q) remaining resolutely centred on the origin. Although the very large autocorrelation times that are 
encountered for (3 > (3sg make the interpretation of the low temperature results a dangerous business we 
can press on deeper into the low temperature region to look at P(q) there and see if there is any evidence 
for the second transition. P(q) is plotted at (3 = 3.5 for the same three lattice sizes in Fig. 9, where we 
can see that there does, indeed, appear to be a second peak developing away from the origin at q = 0.2 
as well as the bump at q = 0.7. 

The qualitative features of the above results are preserved in the Q = 4, 10, 50 state models that 
we also simulated: the histogram of P(q) broadens from a sharp peak at the origin at around (3sg 
and develops a secondary bump or shoulder. There is some evidence of further structure developing 
for even larger (3. From the numerical evidence it would thus appear that the Q = 3 Potts spin glass 
with a ± J distribution of bonds lies above the critical Q value that separates Ising and Potts spin glass 
behaviour, whereas the Bethe lattice calculation in supports the opposite conclusion, although on the 
relatively small lattices simulated here it is difficult to distinguish between features that will be sharp in 
the continuum limit and continuous parts of P(q)- A more detailed analysis of the saddle point equations 
may shed some light on this question analytically. For completeness in Fig. 10 we show P(q) for the Q = 4 
model at (3 = 4.4(> (3sg) to demonstrate the similarity with the Q = 3 results. 

The fits to the maxima of the Binder's cumulant to extract vd are rather poor for all the higher 
state Potts models, but the specific heat curves for Q = 3,4 are similar in form to the Ising model, as 
can be seen for the Q — 4 Potts model in Fig.ll, indicating the continued presence of a cusp rather 
than a divergence. We therefore do not have a reliable value of vd to feed into the finite size scaling 
relations, although the cusp suggests that vd ~ 3, at least for Q — 3,4. Nonetheless, we can still fit to 
find the combinations (3/vd and j/vd at the estimated critical points for the various models. This gives 
(3/vd — 0.28(2), 0.31(5) for the Q — 3,4 state models respectively. The mean field theory suggests that 
q becomes discontinuous at (3sg f° r Q > 4. With the smaller lattices that we have for Q = 10, 50 this 
effect, if it exists, is masked by finite size rounding but attempting to fit (3/vd does give larger values 
and much poorer fits than for Q = 3, 4 such as 0.45(5) for Q = 10. This could be construed as providing 
some evidence for a discontinuity. The data allowed fits to "//vd only for the Q = 3,4 models, where it 
gave 0.36(2), 0.40(7) respectively. 

In summary, all the Potts P(q) are clearly qualitatively different from the Ising model, even though 
the two exponents (3/vd, j/vd we fitted for Q — 3,4 are roughly similar. There is not sufficient data in 
the simulations reported here to reliably distinguish a discontinuous transition for Q > 4, though there 
are certainly indications (the fits to (3/vd) that this is the case. In all the cases, however, the values 
extracted for the spin glass transition temperature by analysis of Binder's cumulant are in agreement 
with the predictions. 

5 Conclusions 

The novel analytical approach offered by using techniques borrowed from matrix models in the manner of 
B to look at spin models on random graphs allows one to rederive results for transition temperatures and 
order parameters that are less transparent in previous replica calculations, or only arrived at by analogy 
with the Bethe lattice. Indeed, results from bifurcation theory enabled us to calculate the Hessians for 
bifurcations from the symmetric solution (which may not be dominant all the way to the bifurcation 
point, as is seen in the k — 3,4.. solutions) in k replica Ising or Potts models for any k and P(J)- The 
number of putative spin glass solutions bifurcating at (3i appeared to behave as in the the infinite range 
model. In essence the investigation of a spin glass on a z— regular random graph is reduced to looking at 
the solutions to the equation 



Further work along these lines is clearly both desirable and possible for the Ising and other Potts models, 
with the nature of replica symmetry breaking in this short range (but still mean field) model being perhaps 




lim k 







(53) 



11 



the most important question. We have made no attempt here to follow the various solution branches 
that bifurcate at /?2, for instance, which would shed light on the low temperature phase, particularly if 
calculations for arbitrary k were still possible. An understanding of the role of the first order solutions 
that appear for k > 2 and their failure to influence the k — * limit is also still missing. We emphasize 
again that the transitions in the thin graph model appear to be identical with those in the Ising replica 
magnet where it is the k = 2 transition that marks the boundary between first order transitions and the 
behaviour seen in the quenched model at k = 0. 

We investigated the models numerically in some detail. Quantitatively, transition temperatures were 
in agreement with those calculated for all the models. The extraction of the critical exponents for 
the Ising model showed clear mean field behaviour, and the p-independence of the spin glass transition 
temperature was also apparent. Qualitatively, the P{q) measured in the various models backed up the 
mean field picture of the phase transition with a continuous distribution for the Ising model and sharper 
features for Q > 3. 

In summary, spin glasses on thin graphs offer a promising arena for the application of ideas from 
matrix models, large-n calculations in field theory and bifurcation theory. The tensor (or near-tensor) 
product of the inverse propagator allows some quite general expressions to be derived for the Hessian 
in the saddle point equations and offers a powerful line of attack on questions such as replica symmetry 
breaking. As a subject for numerical simulations they offer the great advantage of mean field results with 
no infinite range interactions and no boundary problems. 

It is worth remarking in closing that one is not limited to the mean-field theory with the current 
methods. It is possible to "fatten" the graphs analytically by increasing the size of the matrices in the 
saddle point equations N = 2, 3 . . . as has already been done with some success for Ising models coupled 
to two dimensional gravity Jl9[] , which is equivalent to looking at the models on an annealed ensemble of 
planar (fat)graphs. This reduces the fractal dimension of the graphs from infinity in the mean field case 
to more realistic values. Numerical simulations of the ferromagnetic transition |2C| and possible spin glass 
transition Q have already been carried out on a quenched ensemble of such planar graphs and in the spin 
glass case P(q) still presents a mean-field-like appearance. It would be very interesting to say something 
analytically about the nature of the low temperature phase in such a non-mean-field spin glass. 
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Fig. 1 The crossover in Binder's cumulant for the overlap in the Ising spin glass calculated with individ- 
ual disorder averages. Only the 250, 500 and 1000 lattice data is plotted for clarity. The lines arc 
drawn only to guide the eye and are not the best fit curves used to determine the crossing points. 

Fig. 2 The crossover in Binder's cumulant for the overlap calculated with an overall average. Again, 
only the 250, 500 and 1000 lattice data is plotted and the lines are to guide the eye only. 



Fig. 
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The specific heat for the Ising spin glass. 




Fig. 


4 


The spin glass susceptibility Xsg- 




Fig. 
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The linear susceptibility xm- 




Fig. 


6 


P(q) vs q for p = 0, 0.3, 0.5 on graphs of 500 vertices at (3 = 1.2, deep in the 


spin glass phase. 


Fig. 


7 

P 


P(q) vs q for the same p on graphs of 500 vertices at (3 = 0.9, close to (3sg- 
= 0, which deviates slightly, is labelled. 


The histogram for 



Fig. 8 P(q) vs q for the three state Potts model at (3 = 2.2. Data from graphs with 500 and 1000 vertices 
are plotted. 



Fig. 9 P(q) vs q for the three state Potts model at (3 — 3.5. Data from graphs with 500 and 1000 vertices 
are plotted. 

Fig. 10 P(q) vs q for the four state Potts model at (3 = 4.4. Data from graphs with 500 and 1000 vertices 
are plotted. 

Fig. 11 The specific heat for the Q — 4 Potts spin glass. As the transition (determined from the Binder's 
cumulant for q) lies above the specific heat peak the values for the larger lattices are only measured 
on this side. 
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